Dynamics of double-well Bose-Einstein Condensates subject to external Gaussian 

white noise 
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Dynamical properties of the Bose-Einstein condensate in double- well potential subject to Gaussian 
white noise are investigated by numerically solving the time-dependent Gross-Pitaevskii equation. 
The Gaussian white noise is used to describe influence of the random environmental disturbance on 
the double-well condensate. Dynamical evolutions from three different initial states, the Josephson 
oscillation state, the running phase and 7r-mode macroscopic quantum self-trapping states are con- 
sidered. It is shown that the system is rather robust with respect to the weak noise whose strength 
is small and change rate is high. If the evolution time is sufficiently long, the weak noise will finally 
drive the system to evolve from high energy states to low energy states, but in a manner rather 
different from the energy-dissipation effect. In presence of strong noise with either large strength or 
slow change rate, the double-well condensate may exhibit very irregular dynamical behaviors. 

PACS numbers: 03.75.Lm, 03.75.Kk, 03.65.Yz 



I. INTRODUCTION 

The atomic Bose-Einstein condensate (BEC) trapped 
in double-well potentials builds up bosonic Josephson 
junction (BJJ) [l|-|5[- Since it exhibits abundant quan- 
tum properties in comparison to condensate in a single 
trap, the double-well condensate has already been inten- 
sively investigated theoretically in the last few years. As 
a BJJ, the double- well condensate can not only display 
dc, ac Josephson effects and the Shapiro effect. It also 
exhibits an quantum nonlinear effect, named the macro- 
scopic quantum self-trapping (MQST) [j§|. O n the other 
hand, quantum fluctuation is believed to give rise to fas- 
cinating influence on above dynamical behaviors [6Ml2|. 
such as collapse and revival of quantum oscillations (g, [D] , 
disappearance of coherence [&, 9] , and destruction of the 
self-trapped state [10| . 

Experimentally the Josephson tunneling and MQST 
in a single BJJ have been observed in 2005 Since 
then, more progress has been made in studying static, 
thermal and dynamical properties of the double- well con- 
densate. Experimental investigation of thermal induced 
phase fluctuations has been reported [l4| ■ Measurements 
of the ac and dc Josephson effects in BJJ have already 
been realized [l5[ . The BJJ system has also been used to 
perform interference-fringe experiments [16| and to inves- 
tigate the crossover from Josephson dynamics to hydro- 
dynamics [l7| . In above theoretical studies, the double- 
well condensate is mainly treated as an isolated system, 
but actually it is coupled to certain thermal cloud and 
subject to environmental distortions in experiments. Dis- 
sipation and noise effects play important roles in under- 
standing properties of BJJs. 

The dissipative effect has been studied by several 
groups [lUill. It is suggested that the MQST state 
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can be destroyed by energy dissipation [18-20.]. Possi- 
ble decoherence caused by dissipation is also discussed 
[2~il [22[ . It is also shown that dissipation could lead to 
enhancement of coherence under specific conditions |23j . 

Here we consider a kind of noise effect on the double- 
well condensate. Noise can be classified as "internal 
noise" or "external noise" with respect to its origins 
[24| . Internal noise comes from the inside fluctuations of 
the system, including quantum and thermal fluctuations, 
or from the exchange symmetry of identical particles 
[25l [2fjj . On the other hand, "external noise" is brought 
about by fluctuations which are not "self-originating" . It 
is induced by the coupling between the system and its en- 
vironment. For cold atomic systems, external noise may 
originate from the magnetic field, laser beams, or other 
externally applied random driving field. The present 
work will focus on this kind of external noise. 

A number of theoretical works have been devoted to 
the influence of external noise on the double-well con- 
densate [H, [13, Hi]. The noise-induced dephasing [13] 
and phase decoherence [28[ have been predicted. These 
works based on the two-mode approximation mainly dis- 
cussed the phase noise. The phase noise is introduced 
by coupling a stochastic fluctuations either to the tun- 
nel amplitude or to the number-imbalance operator. In 
their treatment, spatial information of the external noise 
is ignored. 

The present paper deals with noise due to fluctuations 
of magnetic field or the optical potential. We shall con- 
sider the time- and space-dependent characteristics of 
the external noise by solving the time-dependent Gross- 
Pitaevskii (GP) equation numerically. In Sect. HH we 
describe a noise model in which the noise is simulated by 
spatially distributed stochastic potentials. The concept 
of Gaussian white noise and the algorithm used to solve 
the time-dependent GP equation for the double- well con- 
densate are briefly described. We present in Sect. Mil the 
obtained results and discuss the influence of noise on the 
dynamical behaviors of system. A brief summary is given 
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in the last section. 



II. THE MODEL AND METHOD 

We consider the BEC confined in a double- well poten- 
tial is composed of atoms of mass m with weakly repul- 
sive interaction. Dynamics of the double- well condensate 
is obtained by solving the time-dependent GP equation 
numerically [29, 30]. Wc formulate the external noise by 
an additional potential of stochastic strength fluctuation 
in the space. 

The external noise is modeled as the Gaussian white 
noise potential V n (z,t), which is space- and time- 
dependent stochastically. So V n (z,t) satisfies the follow- 
ing equations, 



(V n (zi,t))=0, (1) 



and 



{Vn^Vniz^t')) = 2D 5(t - t')6(zi - zj), (2) 

where (} denotes averaging over both the space and the 
time. The Dirac delta function in the correlation formula 
makes sure V n (z,t) is "white" noise, Do is the fluctua- 
tion amplitude of noise potential. Apparently, Dq char- 
acterizes the strength or intensity of noise and it can be 
controlled experimentally. 

The noise correlation function in time at a given po- 
sition is defined as g(t — t') = (V n {z i: t)V n (zi,t'))T — 
(V n {zi,t))T ~ S(t — t') where ()t means that the av- 
erage is just taken over time. Under this definition, 
(V n (z 1 t))x must be zero at all given positions. In numer- 
ical calculations, V n (z,t) is produced as a time-ordered 
series, V n {z,ti), so the time average is performed as 



M 



M is actually a finite num- 



ber, and MAt is the total time interval with At = ti+i—ti 
being the step interval of the noise potential series. In 
our calculations, M = 100 is already enough to makes 
sure (V n (z,t))T = 0. Therefore At can be used to eval- 
uate the velocity of noise. Faster noise corresponds to 
smaller At. 

The general three dimensional time-dependent GP 
equation provides an exact and fundamental description 
for our research, which can be formulated as 



. W(r;t) 

1 dt 



-—V 2 + V(r)+gN\^(r;t)\ 
Zm 



i>(v;t), 
. . (3) 

where ip(r; t) is the macroscopic wave function at position 
r and time t, g — ATrh 2 a/m is the nonlinear interaction 
with a being the s-wave scattering length. The external 
potential V(r;t) consists of the double well Vd w {^) and 
the noise potential V n (r;t). Vd w (r) = ^-{uj 2 x rl + w 2 r 2 + 
w 2 r 2 ) + T4exp(— r 2 z jq^) fixes the double-well configura- 
tion, where wAi = x, y, z) is the trap frequencies in the 
direction of i [12|,[l3| and 14 is the barrier height. Exper- 
imentally, the Bose-Einstein condensates with repulsive 



interaction in a Quasi-ID symmetric double-well poten- 
tial can be achieved by splitting one cigar-shaped atomic 
cloud into two separated aligned cigars by a laser beam. 

Set the double wells being along the z direction and 
in x and y directions the strong confinement is exerted 
(lu x = uj y 3> lo z ). Thus the original three dimensional 
condensates shall be reduced into a quasi one dimen- 
sional cigar BECs and Eq. Q is reduced into a one- 
dimensional one. By scaling the length and energy as 
l z = (muj z ) and Huj z respectively, the dimensionless 
parameters are simplified as z — r z /l z , t — tuj z /2 and 

15 = \J<^ 2 + Wy/uj z . Therefore, the reduced one dimen- 
sional GP equation can be formulated as 



Or 



dz 2 



Vdw(z) + 2v n (z;t) + gi D \ip(z\T)Y 



where Vd w — z 2 + 2«bexp(— z 2 /g 2 ), with v^w — 
V dw /(huj z ), v n = V n /(huj z ) and q' Q = q Q /l z . We ob- 
tain the dimensionless interaction parameter g±o — 
gNm(3 / (irlzh 2 ) with TV the total atom number. 

The wave function ip(z; r) satisfies the normaliza- 
tion condition dz\'ijj(z; r)| 2 = 1. The fraction of 
atom number in the left and right well are Ul(t) — 



J^dzfcr)! 2 
spectively [30J . 



and Ur(t) 
0l(t) = 



= f Q °° dz\i>(z;r)\ 2 



re- 



ar ct an 



f"^ dzIm[4>(z-T)]p(z;T) 
/"^ dzKc{4>{z-T)]p(z-T) 



arctan 



the phase in the 

/o™ dzlm[i>(z-T)]p(z;T) 



left well while 0r(t) 
in the right well with the den- 



/ °° dzRcl4>(z;T)]p(zjT) 

sity p(z;t) = i/j*(z;t)i/j(z;t). 

4>+ (z) and 4>- (z) represent the initial ground state and 
the first excited state wave functions for the condensates 
in the double-well. Their linear combinations can be de- 
fined as the left (right) mode: iPl,r(z) = t±i£J^t=i21 . 
They satisfy the orthogonal condition J dzip l{z)iI) r{z) — 
0. The trial wave function for obtaining the initial state 
can be chosen as the superposition of iPl(z) and 4>r( z ) 
as in Ref. ^(z ;t) = ^ l {t)4> l {z) + ^ b (t)0«(z), 
where i^l(r)(t) = ^/n L ^ R ){T)e ieL ^ R ^ T \ The population 
imbalance and relative phase at time t are defined as 
An = til — riR and A9 — 9l — &r- A given trial at 
r - is #(z;0) = e iA9 ^ ^j0)^ L (z) + yfcrf$il> R {z). 
An(0) = ni(0) — tir(0) represents the initial population 
imbalance and A6(0) the phase difference of condensates 
in double well. The relative phase can be measured by 
the interfere patterns of releasing the Bose-Einstein con- 
densates from the double-well potential after different 
evolution times [l3j]. Moreover, a technique based on 
stimulated light scattering has been developed to detect 
A8 nondestructively [3l|. 

The reduced time-dependent GP equation can be 
solved using the Split-Step Crank-Nicolson scheme, with 
both the space and time being discretized 32(. In our 
calculation, time step is St = 0.001 and space step is 
8z = 0.01. The initial ground state and the first excited 
state wave functions, 4>+(z) and <f>-(z), can be obtained 
by the imaginary-time propagation. The dynamical evo- 
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FIG. 1: Snapshots of the random Gaussian white noise po- 
tential with the noise strength D — 0.005 at the time r = 1 
(a), t = 2 (b), and r = 10 (c), respectively. 



lution are calculated by solving Eq. (|4]) using real-time 
propagation method. 

The Gaussian white noise is produced numerically by 
the Box-Mueller algorithm (33[, 



v n = y/ -AD ln(a) cos(27r&), 



(5) 



where a and b are two uniformly distributed random 
numbers on an unit internal, and the dimensionless pa- 



rameter D = Dq/(T%uj z Y 



varies occasionally from 



time to time. Figure [T] displays snapshots of the spa- 
tial form of the random potential at different given time. 
Suppose that it changes R times each unit time. Then 
R defines a character parameter with respect to the ve- 
locity of noise and it is called the change rate of noise 
hereinafter. The fastest noise is the one whose step in- 
terval At just amounts to the calculation time step St. It 
means that the noise potential changes once per calcu- 
lation step. Since the time step St = 0.001, the change 
rate of the fastest noise in our study is R = 1000. 



III. RESULTS AND DISCUSSIONS 
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FIG. 2: Evolution trajectories in A6-An phase space of a 
system initially in the Josephson oscillation state (An = 0.9, 
A9 = 0) with the fixed change rate of noise, R = 1000. The 
evolution time is from r = to r = 5000 and the noise 
strength is D = (a), D = 0.005 (b), D = 0.02 (c), and 
D — 0.5 (d), respectively. The black horizontal band in panel 
d consists of many almost vertical lines and it denotes the ir- 
regularly energy-fluctuating states in the final evolution stage. 
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FIG. 3: Evolution trajectories in A9-An phase space of a 
system initially in the running phase MQST state (An = 0.9, 
A9 = 0.7n) with the fixed change rate of noise, R — 1000. 
The evolution time is from r = to r = 5000 and the noise 
strength is D = (a), D = 0.005 (b), D = 0.02 (c), and 
D — 0.5 (d), respectively. The black horizontal band in panel 
d describes the same evolution stage as in Fig. [2jd). 



The phase-space diagram of double-well BECs with- 
out noise has been studied based on the time-dependent 
GP equation [2(| and the obtained results are consis- 
tent quantitatively with the two- mode model results §■ 
Three typical dynamical regimes can be present, (i) The 
Josephson oscillation regime consists of a cluster of close 
orbits circling the lowest energy point (An = 0, AO = 0) 
or (An = 0, A9 = 2kir) where k is an integer, (ii) The 
7r-mode MQST regime consists of a cluster of close or- 
bits circling the highest energy point (An ^ 0, AO = 
(2k ± 1)tt). (iii) The running-phase MQST regime con- 
sists of open orbits which lie between the above two 
regimes. 

To obtain the numerical results of dynamical double- 
well condensate subject to the Gaussian white noise, we 
set the double- well potential in Eq. (Q| to be Vdw(z) = 



z 2 + 10cxp(— z 2 ) for simplicity. The interatomic interac- 
tion gi£>=0.01. Under above choice of relevant param- 
eters, all three dynamical regimes appear in the phase 
diagram. First, we reproduce the phase diagram without 
noise, D = 0, and determine the location of each typical 
dynamical regime. Then we choose one initial state from 
each regime by setting the initial particle imbalance An 
and relative phase AO appropriately. 

We first consider the dynamical evolution subject to 
Gaussian white noise with high change rate of R = 1000. 
Fig. (a) shows the orbits presented by constant energy 
lines for an initial state of Josephson oscillation (An=0.9, 
A0=0) without noise (D = 0). Both the population im- 
balance and the relative phase oscillate around the zero 
point (A6 = 0, An = ) during the evolution. When 
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Gaussian white noise with the strength D = 0.005 is im- 
posed on, as illustrated in Fig. &b), the trajectory gets 
a little "fat", but dynamical properties are not qualita- 
tively different from the system without noise. However, 
the noise effect shows up apparently when D — 0.02, as 
seen in Fig. EJc). The system evolves gradually along 
an inward spiral path which is still smooth at earlier 
stage but then becomes zigzag after a long time. Dur- 
ing this period, the system undergoes energy dissipation. 
Fig. [2Jd) shows that the energy dissipates more quickly 
as the noise becomes extremely stronger (D = 0.5). At 
last, the original Josephson oscillation has been com- 
pletely destroyed and the system enters into a state in 
the black horizontal band along the An = line indi- 
cated in Fig. [5Jd). The particle imbalance An is almost 
zero which means that particles tend to distribute almost 
equally in each well, while the relative phase varies very 
fast, signaling that A9 is no longer a well-defined param- 
eter. It is also worth noting that the energy does not 
decrease monotonically any longer, but becomes irregu- 
larly fluctuating. Therefore, such a state is called the 
irregularly energy-fluctuating state in this paper. And it 
is also the final state of the system after a long enough 
evolution, no matter how weak the noise strength is. 

Fig. [3J and 0] demonstrate situations of the initial state 
being in the running phase MQST regime (An=0.9, 
A<9=0.7tt) and 7r-mode MQST (An=0.9, A0=O.9tt), re- 
spectively. In this case, each MQST regime splits into 
two separate parts which are symmetric with respect to 
the An = line in the phase space. 

Look at the running phase MQST case. When the sys- 
tem is subject to a weak noise, the phase may run for a 
few periods and then move around a energy maximum 
point whose evolution trajectory looks like the 7r-mode 
MQST. If the noise strength turns stronger, the sys- 
tem starts 7r-modc-like oscillating from very beginning, 
as shown in Fig. |3J(b). One interesting phenomenon is 
that the trajectory may pass across the An = line and 
then evolves around the down energy-maximum point 
(An < 0, A9 = 7r). This process will be accelerated if D 
is increased and the system may be in a Josephson-like 
oscillation around the point (An=0, A9—tt), as shown 
in Fig. [31(c). The energy decreases with time in an os- 
cillatory manner. In case of extremely strong noise, the 
system falls down to the irregularly energy-fluctuating 
state after a short period of damping. 

Fig. [4] displays the cases evolving initially from the tt- 
mode MQST state (An=0.9, A6»=0.9tt). It is obvious 
that the 7r-mode MQST state is more sensitive to the 
noise. As shown in Fig. 0|b), with the participation of 
noise as weak as D = 0.005, the previous closed orbit 
has already been becoming very fuzzy. When the noise 
strength is increased to D = 0.02, the system gets into 
running phase regime quickly. The following evolution 
can be regarded as a process starting from a running 
phase MQST state. The basic features should be similar 
to Fig. [3] The actual evolution trajectory depends on the 
concrete initial state. This is the reason why Fig. @Jc) 
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FIG. 4: Evolution trajectories in A8-An phase space of a 
system initially in the 7r-mode MQST state (An = 0.9, A6 = 
0.97r) with the fixed change rate of noise, R = 1000. The 
evolution time is from r = to r = 5000 and the noise 
strength is D = (a), D = 0.005 (b), D = 0.02 (c), and 
D — 0.5 (d), respectively. The black horizontal band in panel 
d describes the same evolution stage as in Fig. [21(d) . 




FIG. 5: Evolutions of an initially running phase MQST state 
(An = 0.9, A9 = 0.7tt) in presence of noise with the noise 
strength D = 0.005 and the change rate R = 10. The trajec- 
tory is plotted for the evolution time from r = to r = 5000 
(a), from r = to r = 1000 (b), from r = 1000 to r = 3000 
(c), and from r = 3000 to r = 5000 (d), respectively. The 
An-r plot (e) and Ad-r plot (f) are shown in the third row. 

andQJd) look quite different from Fig. [31c) and[3Jd). 

According to discussions above, when the noise 
strength is small, the Gaussian white noise may lead to 
an energy dissipation effect. The noise with a larger D 
makes the energy dissipating faster. But there are some 
differences in the dissipative process between the present 
study and the results obtained within phenomenological 
theories [l9|, [2(|. The main difference lies in the evo- 
lutions from MQST states. The system damps directly 
towards the saddle point (An=0, A9—tt), not the energy 
minimum (An=0, A9=0). Then it steps into the irreg- 
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ularly energy-fluctuating state. The phenomenological 
theories suggest that the energy decreases all the time 
until the system approaches the energy minimum point. 
No matter how strong the dissipation is, the system will 
finally damp into the energy minimum state, then the 
evolution stops. The irregularly energy-fluctuating state 
does not appear. 

We have calculated the dynamical behaviors of the 
double-well condensate subject to Gaussian white noise 
with different strength D and change rate R. The noise 
effect is mainly determined by the strength D when the 
change rate R is large enough. Nevertheless, the change 
rate R also plays an important role in understanding the 
dynamical evolution when R is small. 

Figure [5] displays a typical case of the noise with small 
strength (D = 0.005) and slow change rate (R = 10). 
The evolution is from an initial state of running phase 
MQST (An=0.9, A6>=0.7tt). Figure [5Jb) looks like 
Fig.EJb), but the evolution time is only within r = 1000, 
far shorter that r = 5000 as in Fig. [SJb). Figure [SJd) 
shows the trajectory with the time from r = 3000 to 
r = 5000. During this period, the system oscillates in 
a Josephson-like way around the saddle point (An=0, 
A0=7r), similar to the case shown in Fig. |3|c). It seems 
that decreasing the change rate gives rise to similar ef- 
fect as increasing the noise strength. The A(9-r and An-r 
plots shown in Fig. [5je) and (f) provide detailed infor- 
mation of the evolution process. Atoms can be trapped 
almost in one well for a while and then oscillate between 
the two wells. The variation range of A9 is significantly 
enlarged and many extremely sharp peaks appear in the 
A9-t line. These peaks indicate that the relative phase 
can change extremely fast, which may signal that the rel- 
ative phase is no longer well-defined, as discussed above. 
Such situation occurs preferently in the large An region, 
as indicated by the almost vertical lines in Fig.[5ja). This 
is different from the large D case, in which the almost 
vertical lines appear in the An = region first. If the 
system evolves further, more almost vertical lines appear 
and they tend to cover the whole phase space in an irregu- 
lar manner. As a result, it is almost impossible to expect 
in which state the system is at a given time. Meantime, 
the energy of the system does not simply tend to damp- 
ing, but may go up and down occasionally. So the energy 
dissipation effect seems smeared out during this stage. 

IV. CONCLUSION 

In conclusion, we have investigated how an external 
Gaussian white noise affects the dynamics of condensates 



in double- well. Dynamical evolutions from three typi- 
cal dynamical regimes of Josephson oscillation, 7r-mode 
MQST and running phase MQST are discussed. It is 
shown that the system keeps its original dynamical fea- 
ture for a long time when it is subject to a noise with 
weak strength and rapid change rate which we called the 
weak noise. In this case, the noise induced energy dissi- 
pation effect is observed. Energy of the system decreases 
with evolution monotonously in the Josephson oscillation 
regimes, and oscillatorily in the MQST regimes. The 
difference between results of the present study and the 
phenomenological theory is discussed. 

Either increasing the noise strength D, or decreasing 
the noise change rate R may give rise to significant influ- 
ence on dynamics of the system. If R is large enough, the 
influence of D is quite clear. The larger D is, the faster 
the state evolves and damps. Moreover, we also figure 
out that the noise with large strength drives the sys- 
tem finally into the irregularly energy-fluctuating state, 
in which the population imbalance tends to being very 
small and the relative phase is no longer well-defined. 
This phenomenon is different from the dissipation effect 
which makes the condensates evolve from the high en- 
ergy state to the lowest energy state where the particle 
imbalance is zero and the relative phase is 2kir (k is an 
integer). Decreasing R brings about similar effect as in- 
creasing D in the early stage of evolution. After a period 
of time, the noise with slow change rate compels system 
to change irregularly, with the relative phase being de- 
stroyed. But this phenomenon happens preferently in the 
large An region in presence of the noise with slow change 
rate, which is different from the effect caused by the noise 
with large noise strength but fast change rate. 
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